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Abstract. The forward neutron production in the ep collisions at 300 GeV measured by HI and ZEUS 
Collaborations at DESY has been used to estimate the total probability for proton fluctuation into tvk + 
and p-K . The probability found is on the order of 30%. This number is compared with the numbers obtained 
for the probability of quark fluctuation into 7r + from several alternative DIS processes (Gottfried sum rule, 
polarized structure function) and the axial- vector coupling constant, where the pion fluctuation is believed 
to play an important role. 

PACS. 12.38.-t Quantum chromodynamics - 12.39.Jh Nonrelativistic quark model 



1 Introduction 

Measurements of the inclusive deep inelastic scattering 
(DIS) have allowed to extract important information on 
the partonic structure of the nucleon. Parton distributions 
have been determined and scaling violations have been 
tested to a high level of accuracy. The QCD parton model 
has been shown to be very reliable in the presence of the 
hard scale. 

On the other hand, semi-inclusive reactions with elec- 
. — ' tromagnetic probes are less explored. Through the study 
k> , of new observables characterizing these processes we may 
j_j ■ ask more detailed questions about the hadronic structure. 
C3 ' In semi-inclusive reactions the hadronic character of the 
nucleon becomes apparent, soft QCD becomes important. 
The most prominent effects in semi-inclusive DIS are the 
consequence of the strong correlation between the sea- 
quark and the sca-antiquark in the relative state with the 
quantum number of the pion £Q . We will refer to this cor- 
relation, as is the custom in the hadronic physics, as the 
pion fluctuation of the nucleon, p — > Nn and p — ► An. 
A similar fluctuation of the proton to K + A, which is less 
strong than the pion fluctuation because of higher exci- 
tation, is the source of the strange sea in the proton and 
explains the ss asymmetry in the neutrino DIS. Also the 
proton fluctuation to DA C may be the source of intrinsic 
charm in proton 

Two collaborations, HI and ZEUS at HERA, have 
measured the production of forward neutrons in electron- 
proton collisions at a center-of-mass energy of about 300 
GeV |415IBI7ftl5| . These events arc successfully described 
within the models based on the pion exchange mechanism: 
the DIS takes place on the 7r + fluctuation of the proton. 



Neutrons carrying more than 50% of the proton initial en- 
ergy produced in DIS on pions can be well distinguished 
from the neutrons produced in ordinary DIS on the proton 
by their energy and four-momentum distributions. 

The forward neutron production with the aim to ob- 
serve the pion fluctuation has already been studied pre- 
viously in the high-energy pp-reaction at ISR in CERN 
|10j . But the rescattering of proton on the virtual neu- 
tron is so strong that a reliable estimate of the probabil- 



ity of the p 



n fluctuation was not possible. On the 
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other hand, the rescattering of the highly virtual photon 
(Q 2 > 10 GeV 2 ) on neutron is sufficiently small so that 
a reliable probability for the pion fluctuation can be de- 
duced. Thus, the forward neutron production in DIS is 
probably the most direct way to determine the probabil- 
ity of the pion fluctuation quantitatively. The pion fluctu- 
ation is very likely responsible for the flavour asymmetry 
of the quark sea in the nucleon as observed in the viola- 
tion of the Gottfried sum rule. Furthermore, the reduction 
of the amount of the angular momentum carried by the 
quark spins in the nucleon is probably also, to a great 
extent, a consequence of a pion fluctuation, as we will 
discuss later. These two phenomena as well as the value 
of the axial current coupling constant are directly related 
to the pion fluctuation of the quarks, u — > d7r + ,it7r° and 
d — > uir~,dir°. The pion Compton length is larger than 
the nucleon radius and the fluctuating pions overlap and 
interfere. They result in different Fock components of the 
nucleon with one, two and three pions. 

Therefore we felt it interesting to estimate the prob- 
ability of the pion fluctuation in the proton in DIS and 
compare it with the probability of the fluctuation of the 
quark. For the neutrons carrying large fraction of the ini- 
tial proton energy the energy spectra agrees well with the 
pion exchange model predictions. At lower energies contri- 
butions from ordinary DIS are becoming important. The 
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HI and ZEUS measurements were fitted by theoretical 
pion fluxes from which the pion probability in the nu- 
clcon can be deduced. Our present analysis investigates 
the uncertainty of this probability as a consequence of the 
different form of the neutron spectra as a function of the 
coupling constant and the formfactor. The contribution of 
the DIS on pion at low neutron energies was calculated by 
extrapolating the neutron spectrum from energies above 
50% of the proton incoming energy to low energies and 
to full range of four-momentum transfer squared (t) using 
three different pion fluxes which fitted the experimental 
data best. 



2 Estimate of the neutron flux from DIS on 
pion 

To evaluate the probability of pion fluctuation in DIS 
we use the measurements of the leading neutron produc- 
tion at HERA performed by HI and ZEUS Collabora- 
tions. There the leading neutron production was studied 
in different processes: in the semi-inclusive DIS process 
ep — > enX, in the photoproduction of D* and in the di- 
jet production 4,5,6,7,8,9 . The kinematic range of lead- 
ing neutron measurements are restricted by the geometri- 
cal acceptance of the forward neutron calorimeters to the 
range xl ~ 0.2 and 9 n $ 0.8 mrad, where xl is the neutron 
energy relative to the initial proton energy and 6 n is the 
angle of the scattered neutron in the laboratory frame. 

In our analysis we use only the measurements of the 
semi-inclusive channel ep — > enX which have the high- 
est statistical significance and for which one expects the 
least distortion providing the measurement is done by suf- 
ficiently high Q 2 . At low Q 2 the distortion can arise from 
the rescattering of the neutron on the extended hadronic 
photon the size of which increases with decreasing Q 2 . 
The rescattering effects are expected to be small in the 
DIS regime, at sufficiently large photon virtualities Q 2 , as 
is also confirmed by the measurements 0. 

Previous studies have shown that particle exchange 
models describe the leading neutron production data both 
in deep inelastic scattering at HERA and at hadropro- 
duction experiments. In these models the transition 
amplitude for p — > n is dominated by n + exchange 12 
IT3lll4lll5| . The interpretation of this process depends on 
the reference frame. In the rest frame of the proton DIS 
takes place from a pion emitted from the proton. In the 
infinite momentum frame (the rest frame of the photon) 
the photon emites the pion. Such ambiguities are resolved 
using light-front Fock methods fT|)J. Let us first consider 
that the production of leading neutrons in DIS at large xl 
proceeds entirely viap — > n + n + channel. Then the cross- 
section for photon-proton scattering to the final state nX 
takes the form 



the hard photon-pion interaction. The general form of 
pion flux is given by the expression 



^ 7VnX = / I +/ Jx L ,t)-da~<* 



(1) 



where / ff + / p (xl, t) is the pion flux associated with the 
beam proton and da 1 
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denotes the cross-section for 



(2) 

Here is the pion mass and g 2 n /4ir — 13.7 is the pirn 
coupling constant, known from phenomenological analyses 
of low-energy data fjjj . G is a formfactor which accounts 
for off-mass-shell corrections, normalized to be unity at 
the pion pole. For the light-cone approach G becomes 
also dependent on xl- The probability of pion fluctua- 
tion depends thus on xl and t. Our final result, called 
< nir + \p > 2 , corresponds to the probability integrated 
over x L and t [HIT??] . 

Several functional forms of pion flux can be found in 
the literature. These expressions differ in the assumptions 
for the formfactor G and in the assumption of reggeization 
of the pion exchange (e.g. the value of ce(t)). In Figf^we 
show different predictions for the shape of the neutron 
energy distribution. The difference in the normalization 
should not affect our estimate as our analysis depends 
only on the ratios of the phase spaces. The extrapolation 
to the low energies and full t range depends strongly on 
the choice of the flux. 

We quote results using the pion fluxes of |13U14II15) 
w hich are also chosen by the experiments ^OOEIEHHIEI • 

In ZEUS analysis [7] the fraction of leading neutron 
events was measured in the kinematic range Q 2 > 4 GeV 2 , 
and it was found to be 8.0±0.5% for x L > 0.2 or 5.8±0.3% 
for xl > 0.49. A similar result is obtained in the HI 
measurements. The result is however published only in 
the PhD thesis |3U], quoting the fraction of leading neu- 
tron events in DIS to be (7.9± 2 ;£)% for x L > 0.49, \t\ < 
0.5 GeV 2 . In this analysis the inclination of proton beam 
with respect to the nominal value (beam tilt) was not 
properly taken into account. Considering the increase of 
geometrical acceptance for the neutrons by about 20% due 
to the observed proton beam tilt as was done for ref. 9 
gives a corrected value of about (6.3 ± 1.5)%, which is in 
agreement with ZEUS results. Because of the larger errors 
of the value obtained by HI, only the ZEUS result will be 
used in further analysis. 

Using the measured value of (5.8 ± 0.3)% we estimate 
the full contribution of pion exchange to the DIS by ex- 
trapolating the measurement from the measured kinemat- 
ical range to the full phase space. The result of this extrap- 
olation is (9.5±0.5)% using the flux from [T3], (10.3±0.5)% 
for flux [HI and (14.6 ± 0.7)% for flux [TS]. 

In this procedure we assumed that all events with lead- 
ing neutrons in the kinematic range xl > 0.49 originate 
from pion exchange via p — > n + tt + channel. However, 
there are also other processes which can contribute. As 
was estimated in [13] in the kinematic range of the mea- 
surements, the leading neutron production due to p and <i2 
exchanges, to pomeron exchange, or to resonance decays, 
is about 10% of that due to pion exchange. The proton 
transition to A + tt may also contribute to the production 
of leading neutrons. Measurements of p — > n andp — > A ++ 
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reactions at Fermilab ;21 indicate that contribution of A 
channel to forward neutron production is less than 6%. 
This observation agrees with theoretical estimate of A + n 
contribution to be one third to two third of the N + tt 
probability [22] , if we take into account all channels which 
can contribute to forward neutron production and the ac- 
ceptance limitation due to larger angular spread and lower 
energies of neutrons produced from A decay. 

Another effect which one has to take into account is the 
absorption or rescattering. It was pointed out by D'Alesio 
and Pirner JH] that also for the forward neutron produc- 
tion in 7p — > nX reaction even for Q 2 > 10 GeV 2 the 
rescattering of the virtual gamma is not negligible. For 
the xl > 0.5 the average reduction of the cross-section 
(e.g. increase of probability) is about 15%. 

Taking these effects into account the final values for 
the corrected fraction of the DIS on the pion compared to 
the total DIS cross-section is (9.9 ±0.5)% for the flux from 
P3], (10.8 ±0.5)% for the flux from [TJ| and (15.2 ±0.7)% 
for the flux from |15| . The big discrepancy between the 
results calculated using fluxes from [rMT%] and flux from 
[To] is the consequence of the fact that in the flux [T5] 
a constant formfactor was assumed. The discrepancy of 
30% between the fluxes comes from the contribution to 
the extrapolated cross-section from the xl < 0.49 region. 
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Fig. 1. The pion flux as function of xl for different shapes 
of the formfactor integrated over t in the range determined by 
the angular acceptance of the forward neutron calorimeter. 
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Fig. 2. The HI measurement of Ff ,JV /F 7r plotted as a function 
of (3 = x-k for fixed values of Q 2 , where r n is the integrated 
pion flux factor. F% /JV can be interpreted as being the pion 
structure function F% . The data are compared to the GRV-LO 
|23| parameterization of F2 of the proton scaled by 2/3. 
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Fig. 3. The ZEUS measurement of Fg in the range 0.64 < 
xl < 0.82 as a function of x-n for fixed Q 2 . The data are com- 
pared to the GRV 1251 and SMRS I2GI parameterizations of 



3 Probability of pion fluctuation observed in 
DIS 

To estimate the full probability of the p — > 7r + n fluctua- 
tion, one needs to know the relative cross-section of DIS 
on pion compared to DIS on proton. For xl ~ 1 the pion 
is close to mass shell and one expects the ratio of pion 
cross-section to total cross-section to be 2/3. The ratio of 
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F£, measured in HI in the range 0.7 < xl < 0.9, to F% is 
consistent with the ratio of 2/3 (21123] (FigEJ. In this mea- 
surement the flux from |13| was used for normalization. In 
the ZEUS analysis |Hj two different normalizations of flux 
factor were considered. The first method gives ratio of 
to F 2 equal to 0.361. Here the normalization was done 
using the effective flux which was taken from hadronic 
interactions |15| . This result contradicts the expectation 
that for xl = 1 the ratio is 2/3. The second method used 
by ZEUS assumes the ratio of 7*7r to -f*p cross-sections to 
be 2/3 at xl = f - This result is in good agreement with 
the GRV parameterization [2H] of (Fig©. The differ- 
ence in the center-of-mass energies of 7*7r and 7*75 is also 
taken into account as 

a ^ {s{ }~ XL)) = f (1 - xi)**™- 1 = 0.6. (3) 

The probability of the 7r + fluctuation of proton is then 
(16.5 ±0.8)% for the flux from Q3], (18.0 ±0.9)% for [Hj 
and (25.3±1.2)% for ^j- Using the isospin symmetry the 
fluctuation of p — > p + ir° is half of these numbers. The full 
probability of the pion fluctuation of the proton is then 
(27.0 ± 1.2)% for the flux from [Ej, (24.8 ± 1.3)% for [Hj 
and (38.0 ± 1.8)% for Q2|. 

We note that the errors are entirely experimental er- 
rors. The systematic errors of our evaluation are reflected 
in the spread of the numbers for probabilities assuming 
different forms of pion fluxes. 



4 Discussion 

The pion fluctuation of the nucleon in the quark model has 
been so far investigated in inclusive deep inelastic scatter- 
ing: Gottfried sum rule, spin structure functions. From 
these experiments the parameter which defines the proba- 
blity of quark fluctuation into pion, a — (<i7r + |u} 2 , can be 
extracted. 

The importance of the pion fluctuation in DIS was first 
realized in the violation of Gottfried sum rule (GSR) |271 
I28| . The GSR expressed in parton distribution functions 
separated in valence quarks («„, d v ) and sea quarks (u, d) 
reads 

1 r 1 

/ G (0,1;Q 2 ) = -/ dx[u v {x,Q 2 )-d v (x,Q 2 )} 
Jo 

-^j\x[u(x,Q 2 )-d(x,Q 2 )}. (4) 

For the symmetric sea 

/ G (0,l,Q 2 ) = i (5) 

is obtained by noting that the first term in JIJ becomes 
simply the difference between the number of up and down 
valence quarks in a proton, times 1/3. 



The second integral of is straitforward to evaluate, 
if one assumes that the entire sea asymmetry comes from 
the pion fluctuation 

/ dx[u{x,Q 2 )-d{x,Q 2 )} =a{Q 2 ). (6) 
Jo 

and 

la = i - fa. (7) 

The first published results of the GSR by the New 
Muon Collaboration [21] quotes 

I G (0,l;Ql = 4 GeV 2 ) = 0.240 ±0.016 (8) 

This number has been afterwards corrected for the shad- 
owing in the deuteron measurement |29| resulting to be 
I G = 0.216 ± 0.024(ea;p.) ± 0.009(shad.). From this num- 
ber using the probability for the pion fluctuation a = 
0.18 ± 0.05. The fluctuation of the nucleon into A + tt also 
contributes to the correction to GSR with the same sign 
as Ntt fluctuation. 

The model used by Eichten et al. [21] to calculate the 
violation of the GSR we can consider as the lowest order 
of the flavour SU(2)-chiral model of the nucleon. Within 
this model Pirner |30j gives an expression for the integral 
of the spin structure function as a function of a as well as 
the expression for the axial- vector coupling constant qa- 

The integral of the spin structure function I p can be 
written as 

/ P = ^(l-2a). (9) 

I p was measured by EMC |3T] as I r = 0.126±0.010(stat.)± 
0.015(syst.), from which the number a can be obtained as 
a = 0.273 ±0.030. 

The next best experimentally measured spin structure 
function is for the deuterium. It was measured by E143 
[32] as I d = 0.042 ± 0.005. The number a is obtained from 

/u = ^(l-3o) (10) 

as 0.283 ± 0.040. 

We further comment on an alternative estimate of Ip 
and Id- In order to avoid the explicite writing-down of 
the baryon wave functions one assumes the validity of the 
flavour SU(3) for the hyperon octet. In the SU(3) model 
the semi-leptonic weak decays can be well fitted by two 
parameters |33j . the irreducible matrix elements of the 
SU(3) representation, F and D. The flavour SU(3) model 
is liked by theorists as it avoids writing explicitly the wave 
functions, but assures their proper antisymmetrisation. It 
may be of some disadvantage when applied to the spin 
structure functions. The axial-vector coupling constant 
expressed in the model is <?a = F + D, and is fitted to 
the experimentally measured value gA = 1.267 ±0.003. In 
the SU(3) model the spin structure function depends also 
on the ratio F/D. This ratio is determined by optimiz- 
ing the fit of the semi-leptonic decays of the hyperons, all 
of them having AS = 1, and takes into account possible 
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kaon fluctuation of the proton. The SU(3) prediction for 
the integrals Ip and Id using the best known values of F 
and D from Cabibbo et al. [22| corresponds in the SU(2)- 
chiral model to a ss 0.19, supporting our assumption that 
the main contribution to the spin reduction comes from 
the pion fluctuation. 

The axial current coupling constant may give an inde- 
pendent information on the probability of the pion fluctu- 
ation in the proton ground state. It was early realized that 
the experimental value of the axial- vector current coupling 
constant of the weak decays gA can be reproduced if the 
pion coupling to the quarks is also taken into account in 
the axial-vector current (partially conserved axial-vector 
current) 3 1 21 ■ However, the estimate of the pion con- 
tribution is srongly model dependent. In the SU(2) chiral 
model |20I 

9A = {l-a)\. (11) 

From the measured value of §a — 1.267 ± 0.003 and using 
eg [HI we obtain a = 0.240 ± 0.002. 

We want to compare our evaluation of the probability 
of pion fluctuation in the proton (Pao) with values of a 
obtained from previously discussed inclusive processes. In 
our evaluation we used data which explicitly require lead- 
ing neutrons in the final state. The relation between Pn-k 
and a is not obvious, a is the probability of pion fluctua- 
tion on the quark level, and Pn-k is the probablity of pion 
fluctuation on the nucleon level. To connect these quan- 
tities one would need the knowledge of the microscopic 
wave function of the nucleon. 

In Fig 0| our values of the probability of the proton 
fluctuation into n + n + (which is 2/3P/v 7r ) are shown to- 
gether with the values of a, the it-quark fluctuation in 
d + ir + , from different processes. In view of the very crude 
approximation of pion fluctuation being the ground state 
and ignoring the Q 2 dependences, the agreement between 
the different values should be considered as reasonably 
good. The observed similarity of the values of a and Pn-k 
demonstrated clearly that a is not an additive quantity 
PoUTU] . the pion cloud from different quarks overlap and 
interfere. It is also another demonstration of the impor- 
tance of soft interactions in semi-inclusive processes. 
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